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We introduce and theoretically study the plasmon-photogalvanic effect in the planar noncen¬ 
trosymmetric plasmonic crystal containing a homogeneous two-dimensional electron system gated 
by a periodic metal grating with an asymmetric unit cell. The plasmon-photogalvanic DC current 
arises due to the two-dimensional electron drag by the noncentrosymmetric plasmon modes excited 
under normal incidence of terahertz radiation. We show that the collective plasmon modes of the 
planar plasmonic crystal become strongly noncentrosymmetric in the weak coupling regime of their 
anticrossing. Large plasmon wavevector (which is typically by two-three orders of magnitude greater 
than the terahertz photon wavevector) along with strong near-field enhancement at the plasmon res¬ 
onance make the plasmonic drag a much stronger effect compared to the photon drag observed in 
conventional two-dimensional electron systems. 


I. INTRODUCTION 

In natural crystals lacking the center of symmetry, a 
DC electric current can be generated under the action 
of optical illumination. This phenomenon is commonly 
referred to as the photogalvanic effect.l^l Physically, the 
DC photogalvanic signal appears as a result of nonlinear 
mixing of the optical fields with frequencies w and —w 
in a noncentrosymmetric nonlinear crystal. Generally, 
the absence of inversion symmetry in a natural material 
yields the asymmetry in the electronic response while the 
optical field may be symmetric or even homogeneous .1^ 
Artificial noncentrosymmetric microperiodic structures 
with a two-dimensional (2D) electron system can also 
exhibit photogalvanic DC responses due to the asymme¬ 
try of the electronic response.l^f^ Under oblique incidence 
of the electromagnetic wave, the photocurrent can be ex¬ 
cited even in a centrosymmetric homogeneous 2D elec¬ 
tron syste m, including graphene, due to the photon drag 
effect.“i^In this case, the asymmetry of the incident elec¬ 
tric held in the plane of the 2D electron system lifts the 
symmetry constrains for exciting a non-zero photocur¬ 
rent. 

It is also known that a photocurrent can be induced 
on the surface of a centrosymmetric metal material if 
the surface of metal is patterned in the form of diffrac¬ 
tion grating with an asymmetric unit cell.^ The metal 
grating couples light to the surface plasmon-polaritons 
on the metal surface since the reciprocal wavevectors 
of the periodic grating ensure the momentum conserva¬ 
tion between the incident light and the surface plasmon- 
polaritons. In case of the diffraction grating with a bro¬ 
ken symmetry, E^HUl the excited plasmon-polaritons are 
predominantly unidirectional and can effectively drag the 
free carriers which results in a net photocurrent 

While the surface plasmon-polaritons in metals are ex¬ 


cited in visible and near-infrared frequency ranges, the 
frequencies of the^asmons in 2D electron systems,^ in¬ 
cluding graphene,!^ fall into the terahertz (THz) range. 
In contrast to the surface plasmon-polaritons in metals, 
the wavelength of the plasmons in the 2D electron sys¬ 
tems is much shorter (by two and even three orders of 
magnitude) than the wavelength of electromagnetic ra¬ 
diation of the same frequency. Therefore, the THz radi¬ 
ation can be coupled to the 2D plasmons by using the 
grating coupler with the (sub)micron period.^ ^ * ^^ * ^° l The 
two-dimensional electron system with a periodic grating 
can be viewed as a 2D plasmonic crystal (2D-PlCr) since, 
in such a system, the electromagnetic near-held has a 
plasmonic nature and its wavelength is determined by 
the periodicity of the grating. The photocurrent induced 
by THz radiation in 2D-PlCr originates from the pon- 
deromotive (plasmonic) nonlinea rities of the free-electron 
motion in a 2D electron systeirP^l^IIIIl ^ 1^0 00 ^- 

trosymmetry is br oken b y a DC current passing through 
2D electron systerrPSHH or because of geometrical asym¬ 
metry of the grating coupler unit celP^ (see Fig. [^. In 
the 2D-PlCr, the asymmetry of the plasmonic near-held 
in the structure unit cell is essential for an appearance 
of the plasmonic THz photocurrent while the electronic 
response in the 2D electron system can be symmetric 
or even homogeneous. The strength of the plasmonic 
near-held can be considerably enhanced by exciting the 
plasmon resonances in 2D-PlCr. As demonstrated both 
theoreticallji^ and experimentally,^2lH2ll 2D-PlCr with a 
geometrically asymmetric unit cell can exhibit a giant 
photogalvanic response at the THz frequencies. 

As shown in Ref. UHl the plasmonic photoresponse of 
the 2D-PlCr with an asymmetric unit cell is caused by 
two different physical mechanisms (which might be com¬ 
parable in strength), namely, the plasmonic drag and 
plasmonic ratchet effects. Both mechanisms rely on an 
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FIG. 1. Schematic view of the 2D-PlCr with an asymmetric 
unit cell. External electromagnetic wave is incident from the 
top upon the 2D-PlCr at normal direction to its plane with 
the polarization of the electric field across the metal grat¬ 
ing strips. Terahertz photocurrent generates the photovoltage 
t/ph between the open side contacts. 


asymmetric profile of the plasmonic field in the 2D-PlCr 
unit cell. Plasmonic ratchet effect is possible only in 2D- 
PlCr with a spatially periodic 2D electron system and 
can exist also for a centrosymmetric profile of the plas¬ 
monic field (which should however be shifted relative to 
the profile of the spatial modulation of 2D electron sys¬ 
tem though). Although the plasmonic drag effect can 
exist even in a homogeneous 2D electron system, the 
noncentrosymmetry of the plasmonic field in 2D-PlCr is 
required in this case. In general, the asymmetric profile 
of the plasmonic field is a direct consequence of the ge¬ 
ometrically asymmetric unit cell. However, the relation 
between the asymmetry of the 2D-PlCr unit cell and that 
of the plasmonic-field profile is intricate, nontrivial and 
dependent on a particular kind of 2D-PlCr. Although 
the plasmonic field is expected to increase at the plasmon 
resonance frequency, the excitation of the plasmon reso¬ 
nance does not in general guarantee a resonant increase 
in the asymmetry of the plasmonic field. Therefore, it 
has not been fully understood up to now why the asym¬ 
metry of the plasmonic field in 2D-PlCr is so strong to 
cause the ultra-high THz plasmonic photoresponse pre¬ 
dicted in Ref. |25] and experimentally observed in Refs. 

[Mm 

In this paper, we consider the 2D-PlCr formed by a 
homogeneous 2D electron system screened by a periodic 
metal grating with an asymmetric unit cell and focus 
on the noncentrosymmetric differential plasmonic drag 
effect induced in such a system by the THz wave inci¬ 
dent at normal direction upon the 2D-PlCr plane. The 
metal grating is formed by two coplanar interdigitated 
sub-gratings having the same period but different widths 
of the metal strips wi and W 2 - Geometrical asymmetry of 
the unit cell can be varied by shifting one sub-grating rel¬ 
ative to the other in the grating plane. Fig. [l] We show 
that the plasmon resonances with strong asymmetry of 
the near electric field can be excited in the anticrossing 


regime between different plasmon modes of the 2D-PlCr 
which leads to plasmon-photogalvanic current caused by 
the differential plasmon drag. In Sec. H, we obtain the 
formula for calculating the differential plasmon-drag cur¬ 
rent in a homogeneous 2D electron system. Plasmon 
mode excitation in the 2D-PlCr with a symmetric unit 
cell is considered in Sec. HI. Noncentrosymmetric plas¬ 
mon modes and induced plasmon-photogalvanic current 
in 2D-PlCr with an asymmetric unit cell are studied in 
Sec. IV. We summarize the main results in Sec. V. 


II. NONCENTROSYMMETRIC PLASMON 
DRAG. NUMERICAL DETAILS 

We assume that the external electromagnetic THz 
wave with the electric field 

t) = Eq exp[-ia;(t -b z/c)] 

polarized across the metal grating strips (the x- 
direction), is incident normally from vacuum upon the 
grating plane. Here uj is the wave frequency, and c is 
the speed of light. The DC photocurrent density gen¬ 
erated due to the ponderomotive nonlinearities in a ho¬ 
mogeneous 2D electron system by normally incident THz 
wave can be described by the equatiorP^ 

L/2 

-L/2 

where cr{ui) = e^V 2 £)r/[m(l —iwr)] is the dynamic Drude 
conductivity of 2D electron system with — e and m being 
the electron charge (e > 0) and effective mass, N 2 D is 
the sheet electron density, r is the characteristic elec¬ 
tron scattering time, and L = wi + Si + W2 + S2 is 
the grating-gate period. The real-valued electric field in 
the plane of 2D electron system is defined as Ex{x,t) = 
[Ex{x) exp{—iujt) -b c.c.]/2. For open side contacts (see 
Fig. [^, the photocurrent generates the photovoltage 
Uph = WRjx between the side contacts, where R and W 
are the DC resistance and width (across the photocur¬ 
rent) of the 2D electron system, respectively. One can 
rewrite Eq. (1) in the Fourier representation as 

jx = -^^Re[a{uj)]'^qp (^lExJ"^ - , ( 2 ) 

p>0 

where qp = 2p7r/L and 

L/2 

Ex,p=^ J A,(x)e-'«-"dx (p = 0,±l,±2...) (3) 

-L/2 

are the amplitudes of the Fourier harmonics of the 
complex-valued electric field Ex{x) in the 2D electron 
system. It is seen from Eq. (2) that the net D C pho¬ 
tocurrent results from the plasmonic dra^^HHl i^y 
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oppositely travelling Fourier harmonics of the near elec¬ 
tric field in 2D electron system. The non-zero differ¬ 
ential plasmon drag current appears only for a noncen- 
trosymmetric near field in the 2D electron system where 
\Ex,p\^ ^ \Ex-p\^- The sign of the photocurrent is de¬ 
termined solely by the asymmetry of the near electric 
field because all other quantities entering Eq. (2) are 
positive and it reverses under the mirror reflection of the 
near field \Ex,p\'^ —t Therefore, this effect can 

be termed the plasmon-photogalvanic effect because it is 
totally determined by the noncentrosymmetry of the near 
electric field in the 2D electron system. It is worth noting 
that the noncentrosymmetric plasmonic drag is allowed 
even for the normal incidence of THz wave upon the 2D- 
PlCr plane despite the lack of the in-plane component 
of the photon momentum in the incoming THz wave. 
This becomes possible due to an asymmetric profile of 
the plasmonic field with non-zero in-plane component of 
the plasmon momentum. 

We simulate the interaction of the external electromag¬ 
netic wave with the 2D-PlCr shown in Fig. 1 within the 
self-consistent electromagnetic approach based on the in¬ 
tegral equation method described in Ref. 1301 We use 
the material parameters N 2 D = 7.5 x 10^^ cm“^ and 
m = 0.2mo, where mo is the free electron mass, which are 
characteristic for the 2D electron system in AlGaN/GaN 
heterostructures.l^ The dielectric constant of the struc¬ 
ture constituent materials is assumed to be e = 9.5 for 
each layer of the heterostructure. The metal grating is 
separated from the 2D electron system by a barrier layer 
of thickness d = 30 nm. In our simulations, we assume 
the metal grating strips being perfectly conductive and 
infinitely thin. This is a quite justified and commonly 
used assumption at THz (and lower) frequencies where 
the metals are characterized by a high real conductivity 
and hence the electromagnetic field penetrates into the 
metal only by a short skin depth (which is much shorter 
than the electromagnetic wavelength). 

In 2D-PlGr under consideration, there are both 
screened and unscreened parts of 2D electron system. 
Accordingly, the screened and unscreened plasma oscilla¬ 
tions can be, in principle, excited in different parts of the 
2D-PlCr. However, the screened and unscreened plas- 
mons have different dispersion characteristics (typically, 
the frequencies of the unscreened plasmons are much 
higher as compared to the screened plasmon frequencies 
for the same plasmon wavevector) so that the screened 
and unscreened plasmon resonances can be excited in¬ 
dependently in 2D-PlGr. For definiteness, we limit our¬ 
selves to studying the screened plasmon resonances ex¬ 
cited at low THz frequencies. The dispersion relation for 
plasmons in 2D electron system screened by a perfectly 
conductive plane iJ^^ 


LO = 


2'Ke^N2Dk 


me 


( 4 ) 


where k is the plasmon wavevector and e = e[l -I- 
coth(fc(i)]. For kd <C 1, one has coth(fcd) « {kd)~^ and 


the dispersion relation Eq. (3) is reduced to 


oj = k 


2'Ke^N2Dd 


me 


( 5 ) 


III. PLASMON MODES IN THE STRUCTURE 
WITH A SYMMETRIC UNIT CELL 


In order to reveal the characteristic properties of the 
plasmonic spectrum in a 2D-PlCr, we assume quite a long 
2D electron scattering time r = 7x 10“^^ s and start with 
considering the 2D-PlCr with a symmetric unit cell. In 
this case the strips of one sub-grating are exactly centered 
in the middle of intervals between the adjacent strips of 
the other sub-grating, si = S 2 = s. Figure shows the 
plasmon absorption spectrum for the symmetric struc¬ 
ture as a function of the strip width W 2 for fixed values 
of Wi = 2 /xm and s = 1.5 /xm. 
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FIG. 2. (Color online) Plasmon absorption spectrum of the 
2D-PlCr with a symmetric unit cell (si = 82 = 1.5 /xm) de¬ 
pending on the variation of the strip width W 2 for a fixed 
strip width wi = 2 /xm and the electron scattering time 
T = 7 X 10”^^ s. 

In the simplest single-strip approximation we can con¬ 
sider the plasmon modes as being excited only in the area 
under an individual metal strip of the width wi {I = 1, 2). 
Taking the open circuit boundary conditions at the strip 
sides we obtain the plasma resonant frequencies uj{l,ni) 
as given by Eqs. (4) and (5) for the size-quantized plas¬ 
mon wavevector 

ki{ni) = nni/wi (rx; = 1, 2, 3 ...) . 

Therefore, in this lowest-order approximation two types 
of screened plasmon modes with I = 1 and I = 2 are 
excited in 2D-PlGr and the Bloch states are constructed 
by the linear combinations of the single-strip excitations 
with the same I and n/. Under the normal incidence of 
THz radiation only the even modes (with odd indexes ni) 
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are excited. (Note that the profile of the in-plane elec¬ 
tric field of the even plasmon mode in 2D electron sys¬ 
tem across a grating-gate strip is even while that of the 
electron density distribution is odd relative to the strip 
center.) In the next order approximation, a coupling be¬ 
tween nearest neighbors is switched on. The resonant 
coupling between the two different plasmon modes be¬ 
comes important if frequencies, U!i{ni) and uj 2 (n 2 ) with 
ni and 712 of the same parity approach each other. 

The horizontal lines outside the narrow anticrossing re¬ 
gions in Fig. [^represent the frequencies of the plasmon 
modes excited under conductive strips of the width wi 
while the frequencies of the plasmon modes excited under 
conductive strips of the width W 2 are inversely propor¬ 
tional to UI 2 according to the dispersion law (4). The anti¬ 
crossing behavior between the plasmon modes excited un¬ 
der the conductive strips of the two different sub-gratings 
is clearly seen in Fig. As followed from the dispersion 
relation (4), the frequencies of two plasmon modes coin¬ 
cide if the widths wi and W 2 are commensurable, i.e., if 
wi/w 2 = ni/n 2 - At the crossing points, the degeneracy 
of the modes is lifted due to their interaction giving rise 
to the plasmon mode splitting and the formation of two 
coupled modes with the frequencies 


LU± 


U}i{ni) + uj2{n2) ^ 
2 


a;i(ni) - a;2(^2) 
^ 2 


-f t 


2 
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where ti 2 is the effective coupling constant equal to the 
sum of the coupling constants of the plasmon mode in a 
particular strip with the plasmon modes in its left- and 
right-hand side neighbor strips. For one of the coupled 
modes the net dipole moments in the areas under all the 
metal strips have the same polarity whereas those for the 
other coupled mode have opposite polarities under the 
adjacent conductive strips, Fig. [^a). Correspondingly, 
the former mode is bright and the latter is subradiant and 
can become optically inactive (dark) at some frequency 
because of the destructive interference of the dipole mo¬ 
ments oscillating under the strips of the different sub¬ 
gratings. In Fig. the dark coupled mode is seen as a 
blind spot in the low-frequency branch of the split modes. 

The plasmon mode profiles of the resonances shown in 
Fig. [^are symmetric with respect to the vertical plane of 
mirror symmetry of the 2D-PlCr, see Fig. |^a). Because 
the condition Ej. p = E^ _p is satisfied in this case for 
both bright and subradiant plasmon modes, no plasmon- 
photogalvanic current flows in 2D-PLCr according to Eq. 
( 2 ). 

Figure 1^ shows the spectra of the spatial Fourier har¬ 
monics of the plasmonic electric field E^ (a;) near the an¬ 
ticrossing point at frequencies around 1.35 THz for the 
2D-PlCr with a symmetric unit cell. It is seen that dif¬ 
ferent Fourier harmonics resonate stronger in different 
coupled resonances due to different plasmon mode pro¬ 
files in those resonances. The resonance linewidth is de¬ 
fined by the electron scattering in the 2D electron system 
(with the contribution 1 /r into the linewidth of any plas¬ 
mon resonance) and radiative damping of the plasmon 


FIG. 3. (Color online) (a) Normalized in-plane electric field 
in the 2D electron system of the 2D-PlCr with a symmetric 
unit cell at the anticrossing regime between the m = 3 and 
712 = 1 plasmon modes at frequency around 1.35 THz for 
W 2 = 628 nm and r = 7x 10“^^ s. The coupled plasmon mode 
resonating at frequency 1.446 THz is bright while the coupled 
plasmon mode excited at frequency 1.296 THz is subradiant. 
The in-plane electric field in 2D electron system is normalized 
to the electric-field amplitude of the incident THz wave, (b) 
The phases of the electric field of the two different coupled 
plasmon modes as functions of the x-coordinate within the 
2D-PlCr unit cell. The difference between the electric fields 
oscillating at different antinodes of the plasmon modes equals 
TT. The positions of the metal-grating strips are shown by 
thick black bars along the x-coordinate. 


mode. The radiative damping is considerably smaller 
for the subradiant plasmon modes and, accordingly, the 
linewidth of the low-frequency resonance is several times 
narrower than that of the high-frequency (bright) one. 
Intriguing feature, however, is that the intensity of the 
absorption resonance of the subradiant plasmon mode 
can be even higher than that of the bright plasmon mode 
at some frequency close to the blind spot in the subra¬ 
diant plasmon-mode branch. This is a consequence of 
the impedance matching requirements called the equipar- 
tition condition: the strongest absorption occurs when 
the dissipative and radiative contributions to the reso¬ 
nance linewidth are equal.^ The radiative damping of 
the subradiant plasmon mode gradually decreases when 
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FIG. 4. (Color online) The spectra of different spatial Fourier 
harmonics of the in-plane component of the plasmonic electric 
field in the 2D electron system near the anticrossing regime 
at frequencies around 1.35 THz for W 2 = 628 nm and r = 
7 X 10“^^ s in the 2D-PlCr with a symmetric unit cell. The 
amplitudes of the Fourier harmonics are normalized to the 
electric-field amplitude in the incident THz wave. 

approaching the blind spot. Hence, at some point of 
the dispersion curve close to the blind spot, the radia¬ 
tive damping of the subradiant mode becomes equal to 
a relatively small dissipative damping, 1/r, which causes 
the strongest absorption at the subradiant plasmon-mode 
resonance at this point. 


IV. PLASMON MODES IN THE STRUCTURE 
WITH AN ASYMMETRIC GRATING. 

PLASMON-PHOTOGALVANIC CURRENT 

Now we turn to the 2D-PlCr with asymmetric unit 
cells where the conductive strips of width W 2 are shifted 
from the center of the gap between the conductive strips 
of width wi making the inter-strip spacings si and S 2 
different. For certainty, we take si = 1 /rm and S 2 = 
2 /rm. As mentioned earlier, a geometrical asymmetry 
itself does not guarantee strong asymmetry of the near 
field in the structure. For example, the Fourier harmonics 
of the electric near field of a bare metal grating (without 
2D electron system at all) with an asymmetric unit cell 
are highly symmetric.!^ Consequently, the photoexcita¬ 
tion of a plasmon resonance in the 2D electron system 
with an asymmetric unit cell, in general, may not lead 
to an enhancement of the plasmonic field asymmetry if, 
e.g., the plasmon mode is excited only under conductive 
strips of one sub-grating. It follows then that only the 
excitation of coupled plasmon modes resonating under 
the conductive strips of both sub-gratings, together with 
a geometrical asymmetry of the 2D-PlCr unit cell, could 
resonantly enhance the asymmetry of the plasmonic field 
in the 2D-PlCr. 

Figure shows the THz absorption spectrum of the 
2D-PlCr with an asymmetric unit cell depending on vari¬ 


FIG. 5. (Color online) Plasmon absorption spectrum of the 
2D-PlCr with an asymmetric unit cell (si = 1 fj,m, S 2 = 2 fim) 
depending on the strip width W 2 for a fixed strip width wi = 
2 fj,m and the electron scattering time r = 7 x 10“^^ s. 


ation of the strip width W 2 for a fixed width wi. Al¬ 
though this spectrum resembles that in the 2D-PlCr with 
a symmetric unit cell shown in Fig. it is much richer. 
Apart from the even plasmon modes (with odd n;), the 
odd plasmon modes can also be excited under conductive 
strips of either sub-grating due to the asymmetry of the 
2D-PlCr unit cell. The odd modes 0 ^ 1 ( 712 ) for n 2 = 4,6... 
above 1.3 THz are well pronounced. The lowest odd 
mode Wi(2) however, which lies at the frequency 0.9 THz 
(between a;i(l) = 0.6 THz and a;i(3) = 1.3 THz), can be 
barely seen near its crossings with the mode a; 2 (l) at W 2 
= 1.1 /xm and mode a;2(3) at W 2 = 3.3 /xm. Therefore, the 
odd mode uJi (2) is totally driven by the latter even modes 
which act as master oscillators. No anticrossing behavior 
is seen in this case because the odd mode provides only a 
weak feedback on the master (even) plasmon modes. The 
odd mode ^ 2 ( 2 ) is seen only for W 2 < 1.5 /xm. Evidently, 
this occurs because the odd mode is more sensitive to the 
asymmetry of the 2D-PlCr unit cell in shorter plasmonic 
cavities where the electric field antinodes are closer to 
the edges of the conductive strip (for the same reason, 
the higher-order odd modes, having the antinodes close 
to the edges of the conductive strip, are well pronounced 
at high frequencies in the spectra shown in Fig. [^. As 
also seen in Fig. the 7 x 2 = 2 plasmon mode is reso¬ 
nantly excited driven by the rxi = 1 plasmon mode at 
frequency 0.6 THz for W 2 = 3.4 /xm. 

Away from the anticrossing points, the normally inci¬ 
dent THz radiation excites the plasmon modes with the 
electric field confined under the strips of either one or an¬ 
other sub-grating. Therefore, the profiles of the plasmon 
modes remain highly symmetric in the x-direction rela¬ 
tive to the center of the respective strip and, hence, the 
amplitudes of the oppositely travelling Fourier harmonics 
remain equal in this case. In the region of anticrossing 
of two even modes (Fig. [^, the plasmonic field of the 
coupled plasmon modes resides simultaneously under the 
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conductive strips of both sub-gratings. Consequently, the 
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FIG. 6. (Color online) (a) The coupled plasmon mode pro¬ 
files in the 2D-PlCr with an asymmetric unit cell at the an¬ 
ticrossing regime between the ni = 3 and n 2 = 1 plasmon 
modes at frequencies around 1.35 THz for W 2 = 628 nm and 
T = 7 X 10“^^ s. The in-plane electric field in 2D electron 
system is normalized to the electric-field amplitude in the in¬ 
cident THz wave, (b) The phases of the electric fields of the 
two different coupled plasmon modes as functions of the x- 
coordinate within the 2D-PlCr unit cell. The positions of the 
metal-grating strips are shown by thick black bars along the 
2 ;-coordinate. 

total plasmonic field is asymmetric throughout the 2D- 
PlCr unit cell. However, in spite of the asymmetric pro¬ 
file of the total plasmonic field, the Fourier harmonics of 
this field are only slightly asymmetric, \Ex^p\ « \Ex^-p\, 
with the maximum relative difference between the am¬ 
plitudes of the oppositely travelling Fourier harmonics of 
the electric field below 1% and 10% in the bright and 
subradiant modes, respectively. This is because the to¬ 
tal plasmonic field of each coupled mode is a sum of the 
two fields resonantly excited under the two sub-gratings, 
while each of the fields is symmetric in its sub-grating 
and the phase shift between different antinodes of either 
plasmon mode is close to tt (see Fig. [^b)). This allows 
for approximating the total field as a product of a con¬ 
stant phase factor and a real function of x-coordinate in 
each coupled mode. 

A remarkable difference in the Fourier-harmonic am¬ 
plitudes \Ex,p\ and \Ex,-p\ can be obtained for a shorter 


electron scattering time, which corresponds to the weak 
coupling regime. In this regime, the width of each plas¬ 
mon resonance is comparable to (or larger than) the split¬ 
ting of the coupled plasmon modes. As a result, the plas¬ 
monic fields of the two coupled plasmon modes strongly 
overlap in the frequency domain which dramatically in¬ 
creases the asymmetry in the amplitudes of the oppo¬ 
sitely travelling Fourier harmonics. In this case, the dif¬ 
ference — \Ex^-p\ can become comparable with the 
amplitude of the Fourier-harmonic itself in the resonance 
of the subradiant plasmon mode (Fig. [^. This is because 



Frequency (TFIz) 



FIG. 7. (Color online) (a) The spectra of different spatial 
Fourier harmonics of the plasmonic electric field in 2D elec¬ 
tron system near the anticrossing regime between the ni = 3 
and n 2 = 1 plasmon modes at frequencies around 1.35 THz 
for W 2 = 628 nm and r = 5 x 10“^^ s in the 2D-PlCr with 
an asymmetric unit cell. The amplitudes of the Fourier har¬ 
monics are normalized to the electric-field amplitude in the 
incident THz wave, (b) The phases of the electric field of the 
two different coupled plasmon modes as functions of the x- 
coordinate within the 2D-PlCr unit cell. The positions of the 
metal-grating strips are shown by thick black bars along the 
x-coordinate. 

the phase of the resonating coupled plasmon mode is dis¬ 
turbed due to its strong coupling with the off-resonance 
plasmon mode. It is seen from Fig. that the elec¬ 
tric field at different antinodes of the subradiant plas¬ 
mon mode oscillates strongly out of phase. The phase 
is predominantly disturbed in the subradiant plasmon 
mode due to large difference in the strength of the dipole 
moments of the subradiant and bright plasmon modes. 
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Enhanced difference in the amplitudes of the oppositely 
travelling Fourier harmonics at the coupled plasmon res¬ 
onances causes corresponding peaks in the photocurrent, 
Eq. (2), due to the differential plasmon drag (Fig. [^. It 
is worth noting that, despite the asymmetry in the am¬ 
plitudes of the oppositely travelling Fourier harmonics of 
the plasmon field is smaller for longer electron scattering 
time, the plasmon-photogalvanic current is much greater 
due to higher quality factor of the plasmon resonance in 
this case. 



FIG. 8. (Color online) Responsivity peaks of the plasmon- 
photogalvanic current at the coupled resonances of the ni = 3 
and n 2 = 1 plasmon modes at frequencies around 1.35 THz 
(w 2 = 628 nm) for two different values of the electron scat¬ 
tering time. The responsivity is defined as the ratio between 
the photocurrent density and the energy flux in the incident 
THz wave. 


The asymmetry of the plasmonic field in the 2D-PlCr 
is even more enhanced in the frequency region of the 
resonant mixing of the even and odd plasmon modes ex¬ 
cited under the conductive strips of different sub-gratings 
(Fig. I^a)). In this case, the difference in the amplitudes 
of the oppositely travelling Fourier harmonics is remark¬ 
able even for the longer electron scattering time because 
the electric fields of the two coupled modes strongly over¬ 
lap in the frequency domain due to a small frequency 
splitting of these modes (Fig. [^b)). Correspondingly, 
the plasmon-photogalvanic current due to the differen¬ 
tial plasmon drag becomes greater (Fig. 101. 

At the resonant frequencies, the electric field \Ex{x)\ 
in the 2D electron system is greater than the amplitude 
of the electric field in the incident THz wave roughly by 
the plasmon resonance quality factor which is typically 
20-30 for r = 7 X 10“^^ s and 2-5 for r = 5 x 10“^^ s 
at the frequency ~ 1 THz. Due to this fact, the 
plasmon-photogalvanic current in 2D electron system, 
Eq. (2), strongly enhances at the frequencies of the 
noncentrosymmetric plasmon resonances. Another rea¬ 
son for the strong plasmon-photogalvanic current is a 
large value of the plasmon wavevector which is from 
two to three orders of magnitude higher than the pho¬ 
ton wavevector in THz frequency range. As a result, the 
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FIG. 9. (Golor online) (a) The spectra of different spatial 
Fourier harmonics of the in-plane component of the plasmonic 
electric field in 2D electron system for the resonant excitation 
of the ni = 2 plasmon mode driven by the 712 = 1 plasmon 
mode at frequencies around 0.9 THz for W 2 = 1.1 pm and 
the electron scattering time r = 7 x lO”'^^ s in the 2D-PlCr 
with an asymmetric unit cell, (b) The plasmon mode pro¬ 
file at the resonant excitation of the ni = 2 plasmon mode 
driven by the 712 = 1 plasmon mode. The amplitudes of the 
Fourier harmonics and the in-plane electric field in 2D elec¬ 
tron system are normalized to the electric-field amplitude in 
the incident THz wave. The positions of the metal-grating 
strips are shown by thick black bars along the ^-coordinate 
in panel (b). 


plasmon-photogalvanic current exceeds the photon drag 
current in 2D electron system by several orders of mag¬ 
nitude. 


V. CONCLUSION 

In conclusion, we have shown that the plasmon- 
photogalvanic current can be excited in 2D-PlCr with 
an asymmetric unit cell via the differential plasmon drag 
of free carriers in 2D electron system by noncentrosym¬ 
metric plasmonic field. A strong asymmetry of the plas¬ 
mon field arises due to the resonant coupling of differ¬ 
ent plasmon modes in the 2D-PlCr. Although we lim¬ 
ited ourselves to the analysis of the coupling and mixing 




























FIG. 10. (Color online) Responsivity peaks of the plasmon- 
photogalvanic current at the resonant excitation of the ni = 2 
plasmon mode driven by the n 2 = 1 plasmon mode at frequen¬ 
cies around 0.9 THz {w 2 = 1.1 /rm) for two different values of 
the electron scattering time. 


between the screened plasmon modes, our results have 
a more general value because they provide a symme¬ 
try criteria which can also be applicable for building a 
strongly asymmetric plasmonic field in a polymodal 2D- 
PlCr supporting both screened and unscreened plasmon 
modes and incorporating a spatially periodic 2D electron 
plasma. Resonant excitation of the plasmonic field with 
a strong asymmetry gives a key to understanding the ori¬ 
gin of ultra-strong THz plasmon-photogalvanic effects in 
2D-PlCr. 
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